Methods to calculate scalar two-loop vertex diagrams 



oo 

On 



Oh 

O 
43 



J. FLEISCHER 
Fakultat fiir Physik, Universitat Bielefeld 

D-33615 Bielefeld, Germany 
E-mail : fleischer @physik . uni-bielefeld . de . 

M. TENTYUKOVQ 
Joint Institute for Nuclear Research, 
141980 Dubna, Moscow Region, Russian Federation. 
E-mail : tent ukov@ thsunl . j inr . dubna. su 



Abstract 
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^\ ' We present a review of the Bielefeld-Dubna activities on multiloop calculations. 



In the first contribution of the above authors in these proceedings, we have introduced our system 
for the automation of evaluation of Feynman diagrams, called DIANA (Diagram ANAlyser). In this 
P-h ■ contribution methods for the evaluation of scalar two-loop integrals will be discussed. 



1 Expansion of three-point functions in terms of external mo- 
menta squared 

04 [ Taylor series expansions in terms of one external momentum squared, q 2 say, were considered in B, 

Pade approximants were introduced in and in Ref. Q it was demonstrated that this approach can 
be used to calculate Feynman diagrams on their cut by analytic continuation. The Taylor coefficients 
are expressed in terms of "bubble diagrams" , i.e. diagrams with external momenta equal zero, which 
makes their evaluation relatively easy. In the case under consideration we have two independent external 
momenta in d — 4 — 2e dimensions. The general expansion of (any loop) scalar 3-point function with its 
i' momentum space representation C(j>\,p2) can be written as 



C(pi,pa)= E a lmn{pi) l (p 2 2) m (Pm) n (1) 
l.m.n— 

X 

' where the coefficients ai mn are to be determined from the given diagram. 

For many applications it suffices to confine to the case p\ = p\ = 0, which is e.g. physically realized 
in the case of the Higgs decay into two photons (H —>■ 77) with p\ and P2 the momenta of the photons. 
In this case only the coefficients a^ n are needed. 

In the two-loop case we consider the scalar integral {k% = k\ — &2, see also Fig. [|) 

C(mi, • • ■ ,m 6 ;p 1 ,p 2 ) = 

(2) 

1 r d i k 1 d i k 2 

(*7T 2 ) 2 J ((fei+p 1 ) 2 -m?)((fei+P2) 2 -m^)((fc2+Pl) 2 -n^)((fe t +P2) 2 -mI)(fc^-r^)(fe^-r^)- 

k t in line 4 (with mass 777,4) depends on the topology: for the planar diagram we have k t = k 2 while for 
the non-planar we have k t = k 3 . 

With obvious abbreviations for the scalar propagators: c; the i th scalar propagator of (U) with p\ = 
P2 = 0, we can quite generally write for the 77"* Taylor coefficient H|: 



(i7T 2 ) 2 a 00n = —— \ d^k^k^ . (3) 

n+1 J a c 2 c 3 c 4 c 5 c 6 



* Supported by Bundesministerium fiir Forschung und Technologie under PH/05-7BI92P 9. 
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Figure 1: Planar and non-planar scalar vertex diagrams and their kinematics 



The numerator F„ can be written as 



For the planar diagram (k t = k 2 ) we have 



-(n-v) „-v ~-( n - v ') „-v' 



E 

i/'=0 



4U*i>k) = £ aZ(k 2 1 ) n - il/+l/ ' )+ ^k 2 2 r(k 1 k 2 y+»'-^. 

fi— max(0,f — n) 



(4) 



(5) 



The coefficients o™^, are mass independent and have been calculated with FORM up to order e 2 (d 
4 — 2s) and stored for the first 30 Taylor coefficients, i.e. they are given in terms of rational numbers. 
Finally all remaining integrals can be reduced by recursion J5| to bubble-integrals of the type 



d d k l d d k 2 



(kf - m 2 ) a (k 2 - mlf^k! - k 2 ) 2 - m§)T ' 



(6) 



or to factorizing one-loop integrals. The genuine two-loop bubble integrals are reduced by means of recur- 
rence relations to Vm({m}). This can be done numerically for the arbitrary mass case or also analytically 
for special cases like e.g. only one non-zero mass. For details see To perform the recursion nu- 

merically, it is important to use the multiple precision FORTRAN by D.Bailey (||) since tremendeous 
cancellations occur in this case. 

We have presented one approach for the calculation of the Taylor expansion of Feynman diagrams in 
some detail, others were worked out in Rcfs. The latter one is particularly suited for programming 

in terms of a formulae manipulating language like FORM. 



2 The method of analytic continuation 

We assume the Taylor expansion of a scalar diagram to be given in the form C(pi,p2, • • •) = J2m=o a my m = 
f(y) and the function on the r.h.s. has a cut for y > yo- If the diagram has a threshold at 4m 2 it is 
suggestive to introduce y = q 2 /4m 2 with q 2 = (pi — p 2 ) 2 as adequate variable with y n = 1. 

The method of evaluation of the original series consists in a first step in a conformal mapping of the 
cut plane into the unit circle and secondly the reexpansion of the function under consideration into a 
power series w.r.t. the new conformal variable. A variable often used is 

c = i-^E2i. ( 7) 

1 + ^1 - y/va 

By this conformal transformation, the y-plane, cut from y to +00, is mapped into the unit circle (see 
Fig.|^) and the cut itself is mapped on its boundary, the upper semicircle corresponding to the upper side 
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Figure 2: Conformal mapping of the y = g 2 /4m^-complex plane into the w-plane. 



of the cut. The origin goes into the point lu = 0. 

After conformal transformation it is suggestive to improve the convergence of the new series w.r.t. to 
by applying the Pade method (^] , |l0[ . A convenient technique for the evaluation of Pade approximations 
is the £-algorithm of [^) which allows one to evaluate the Pade approximants recursively. 

Many interesting examples for the efficiency of this method have been worked out in the meantime 
and we refer to the literatur for massive diagrams to p2| , |l3[ e.g. 

3 Two-loop vertex diagrams with zero thresholds 

Concerning the vertex diagrams, there are many different topologies contributing to a 3-point function 
in the SM. For our purpose of demonstrating the method, we confine ourselves to the planar case shown 
in Fig. ||a. Fig. [3b,d presents infrared divergent diagrams. 
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Figure 3: Planar diagrams with zero thresholds, a) — generic; b) — IR finite diagram; c),d) — IR divergent 
diagrams. 



As a typical (and very complicated case) we consider here Fig. [|b with equal non-zero masses ]l4| . 
There are two massless cuts so that we shall have the double logarithm in the expansion. After summing 
up four contributions we see that the double and single poles in e cancel as well as the scale parameter 
(i, with the result (x = q 2 /M 2 ) 

oo 2 



n=0 j=0 



M 4 



{fo(x) + fx(x) ln(-z) + f 2 (x) \n 2 (-x)} , 



(8) 



where the fj n are now given in terms of rational numbers and C(2)- /^(^O can be summed analytically, 



3 



Table 1: Results for timelike q 2 for diagram ^d. 



q A /M' A 


[12/12] 




[15/15] 






Rc 


Im 


Rc 


Im 


0.05 


+2.948516245 


20.938528 


+2.948516245 


20.938528 


0.1 


-1.108116127 


16.04132127 


-1.108116127 


16.04132127 


0.5 


-4.820692281 


5.066080015 


-4.820692281 


5.066080015 


1.0 


-3.890154 


1.67549787 


-3.890156 


1.67549788 


1.5 


-2.904588 


0.42979 


-2.904581 


0.429778 


2.0 


-2.18294 


-0.06976 


-2.182981 


-0.069728 


10.0 


-0.191 


-0.208 


-0.194 


-0.215 



yielding 

f 2 (x) = (ln\l + x\ -i7rd{l + x)) 2 /x 2 . (9) 

Thus, we have to Pade approximate /o and /i only. Close to the second threshold at q 2 — M 2 the 
convergence is indeed excellent (see Table [l]) Q . It should be noted that for the physical application we 
have in mind, i.e. Z — > bb, this is just the case of interest. It is worthwhile to note the sharp increase for 
low q 2 due to hi 2 (-q 2 /M 2 ). 

The infrared divergent diagrams (Fig. |^c,d) have been successfully considered in Jig ]. 



4 The Differential Equation Method 

We saw that to obtain the expansion of a diagram one has to go through a rather combersome machinery. 
The more coefficients are asked for, the more efforts and machine resourses are required. Thus it is very 
desirable to have analytic expressions for expansion coefficients whenever possible. This can be done 
with the aid of the Differential Equation Method (DEM) |l(J if only one non-zero mass occurs. The 
DEM allows one to get results for massive diagrams by reducing the problem to diagrams with simpler 
structure. 

Let us introduce a graphical notation for the scalar propagators (in euclidean space-time) 

(?T =—-"-— . ^ + m *)« = 

a and m refer to the index and mass of a line. Then one can derive the following recurrence relation for 
a massive triangle fig] 




+Oi2 



+ (ct2 «-> a 3 ) 



Qi- 1 " — 1 ai Qi 

along with some other graphical relations (see details in |]l6|). 

Using this technique we analysed in |]l7f the class of 3-point two-loop massive graphs. As an example 
for the diagram of Fig. ^3 we get 



v y v 



(D-4)J 5c = 2 >^ -2 





— Am' 




- 2m' 




(10) 



In the r.h.s. of (10) the last two terms can be combined, resulting in dJ /dm 2 while for the second 
term we proceed in turn as 
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(D-A) 




(11) 



Thus we are left with simple diagrams (these can be done completely by Feynman parameters) and 
the derivative of the initial diagram w.r.t. to 2 . The solution of the corresponding differential equation in 
terms of a series obtained from an integral representation reads 



J?,c. — — - 



r 2 (l + e) ^ (_ x )"r 2 (n) 



E 



(<7 2 f(m 2 ) r(2n + l) 



1 



1 



-5 - -[Hx) + S 1 (n-1)\ --S 2 (n-1) 



^-S 2 {n - 1) + 4Si(n - l)Si(2n) - C(2) - ln{x)S 1 {n - 1) + i In 2 (a;) 



where 



n , 



A similar formula was obtained for the diagram of Fig.|d (sec 0) 



5 Large mass expansion versus small momentum expansion 

If there are two or more different masses involved the coefficients ai mn in (Q) are not just numbers any 
more but complicated functions of mass ratios. If one of the masses is large, e.g. if the diagram contains 
a top, one can try to perform a large mass (LM, see Jl8| ) expansion rather than a small momentum 
expansion. 

Let us consider the diagram with toi = — m top , m\ — = = and me = myy (see e.g. pj|). 
There are again three subgraphs to be expanded. By direct evaluation we find that there are induced 
poles of the order 1/e 3 in subgraphs while in the sum they cancel which serves as a good check. For this 
particular diagram the result of the LM expansion looks like 




with L = log(TO 2 op /^ 2 ) and b„ s being known functions of q 2 , and [i 2 . 

The considered diagram is evaluated numerically at fj. = m top = 180 , mw = 80 , q 2 = 90 2 . It turns out 
that also here again it is extremely useful to improve the convergence of the LM expansions by applying 
Pade approximants. Up to a normalization factor j^z^— we obtain for the diagram (-9.996,17.9527), 

taking into account terms up to n=10 in (12). The Taylor expansion with 8 coefficients yields a precision 
of appr. 10 decimals the reason for which is that q 2 is far below the threshold. The advantage of the LM 
expansion, however, is that it is easier to program and for many purposes the achieved precision is high 
enough. 
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6 Conclusion 



Any involved calculation in field theory nessesarily consists of two parts: l)automatic generation of 
Feynman diagrams and source codes and 2)techniques of evaluating scalar Feynman diagrams. In our 
contributions we presented both and the methods have a standard that soon important applications can 
be considered. 
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